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ABSTRACT 


The  compartment  heating  system  on  the  U.S.  Coast  Guard’s  Icebreaking  Tug 
(WTGB)  class  cutter  was  studied  to  determine  heat  transfer  performance  characteristics 
of  existing  heat  exchangers  when  used  with  circulating  hot  water  vice  steam. 
Characterizations  such  as  Reynolds  number  vs.  Colburn  j  factor  plots,  convection 
coefficients,  overall  coefficients,  and  Effectiveness-NTU  relations  were  generated.  Initial 
analysis  with  acknowledged  conservative  definitions  of  air  side  convection  coefficients 
determined  that  the  hydronic  system  provided  on  average  seventy  percent  of  the  heat 
transfer  capabilities  available  with  the  steam  system.  Improvements  to  the  hydronic 
system  were  shown  to  increase  heat  exchanger  performance  parameters  by  an  average  of 
ten  percent.  It  was  notable  that  the  added  heat  transfer  available  from  steam  is  not  due  to 
a  property  of  steam  itself  such  as  latent  phase  change  effects,  but  is  due  solely  to  the 
increase  in  entering  tube  side  temperature.  Judging  by  heat  transfer  capabilities  alone, 
with  the  described  conservative  assumptions  on  which  these  results  are  based,  use  of 
currently  installed  heat  exchangers  in  a  hydronic  system  is  a  viable  option. 
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1.  INTRODUCTION 


A.  CUTTER  OPERATIONS 

The  United  States  Coast  Guard  operates  nine  Icebreaking  Tug  (WTGB)  class 
cutters.  Stationed  in  the  northeast  United  States  and  on  the  Great  Lakes,  these  cutters  are 
primarily  used  for  domestic  icebreaking,  but  also  routinely  perform  other  missions  such 
as  search  and  rescue,  pollution  response,  law  enforcement,  and  aids  to  navigation  support. 
The  cutters  are  140  feet  long,  have  a  beam  of  37.5  feet,  displace  662  tons,  and  are 
typically  crewed  by  17  personnel  (see  Appendix  A  for  cutter  illustration).  The  twin 
diesel-electric  propulsion  plant  has  a  cruising  range  of  4,000  miles,  maximum  speed  of 
14.7  knots,  and  provides  sufficient  power  for  the  hull  to  break  through  18  to  20  inches  of 
ice.  A  very  capable  and  versatile  platform,  these  cutters  and  crews  make  vital 
contributions  to  the  overall  service  the  Coast  Guard  provides  to  the  public. 

B.  COMPARTMENT  HEATING  SYSTEM 

1.  Purpose 

The  primary  purpose  of  the  compartment  heating  system  is  to  uphold  the  cutter’s 
overall  mission  readiness,  keeping  on  board  personnel  physically  fit  and  mentally  alert  by 
providing  an  atmosphere  of  suitable  air  for  breathing  under  conditions  that  will  enable  the 
body  to  maintain  a  proper  heat  balance.  A  secondary  purpose  of  the  system  is  provide 
suitable  temperature  and  humidity  conditions  for  the  preservation  of  stores  and 
equipment.  [Ref.  1] 
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2.  Description  of  Present  Configuration 

Compartment  heating  is  accomplished  either  by  duct  heaters  in  the  supply  air 
ventilation  system  or  by  unit  heaters  mounted  in  the  compartments  themselves.  These 
duct  and  unit  heaters  are  either  steam  or  electric.  The  scope  of  this  thesis  includes  study 
of  the  steam  duct  and  unit  heaters  only. 

Heating  of  ventilation  supply  air  is  accomplished  in  two  stages;  first  by  a 
preheater  and  then  by  a  reheater.  Preheaters  are  usually  located  near  supply  air  ventilation 
inlets,  and  heat  incoming  air  sufficiently  to  prevent  condensation  in  ventilation  ducts. 
Reheaters  are  located  in  the  compartments  being  heated,  and  further  heat  the  air  to  a  set 
room  temperature.  Unit  heaters  are  installed  in  compartments  which  require  a  spot  source 
of  heat.  [Ref.  1] 

3.  Deficiencies  of  Present  Configuration 

The  present  steam  system  in  recent  years  has  become  very  maintenance  intensive, 
particularly  with  the  auxiliary  boilers.  The  current  boilers  are  over  15  years  old,  are 
“water  tube”  type,  and  are  difficult  and  expensive  to  procure  spare  parts  for.  With  many 
cutters  operating  in  demanding  harsh  winter  conditions,  reduction  of  heating  equipment 
down  time  is  a  priority.  The  steam  heating  system  equipment  is  also  non-standard 
compared  to  other  ship  classes  in  the  fleet,  leading  to  support  and  maintenance 
difficulties. 
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4. 


Considerations  for  a  Replacement  System 


The  search  for  a  suitable  replacement  for  the  boilers  has  led  to  consideration  of  a 
hydronic  (circulating  water  in  a  closed  piping  system)  compartment  heating  system  to 
replace  the  current  steam  system.  Major  factors  in  choosing  a  replacement  system  are 
ease  of  installation  and  affordability.  A  way  to  ease  the  installation  and  decrease  the  cost 
of  a  new  system  is  to  leave  in  place  and  utilize  as  much  of  the  equipment  from  the  old 
system  as  possible.  There  is  therefore  considerable  motivation  to  determine  the  heat 
transfer  characteristics  of  the  presently  installed  heat  exchangers  when  used  with 
circulating  hot  water. 
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II.  HEAT  EXCHANGER  DESCRIPTIONS 


A.  SYSTEM  OVERVIEW 

The  compartment  heating  system  consists  of  a  steam  generation  plant  and  various 
heat  exchangers.  Included  in  the  steam  generation  plant  are  two  auxiliary  boilers, 
feedwater,  chemical,  and  condensate  tanks,  and  associated  pumps,  controls,  piping,  and 
valves.  The  boilers  are  rated  at  620  poimds  of  saturated  steam  per  hour  at  a  working 
pressure  of  35  psi  with  feedwater  entering  at  210  degrees  F.  A  functional  diagram  of  the 
steam  heating  system  is  shown  in  Figure  1 . 

B.  HEAT  EXCHANGERS 

1.  General 

The  heat  exchangers  installed  on  board  the  WTGB  cutter  class  fall  into  two 
categories;  unit  heaters  and  duct  heaters.  Unit  heaters  consist  of  a  coil,  fan,  motor,  and 
casing  assembly.  They  are  installed  in  machinery  and  work  spaces  throughout  the  cutter. 
Duct  heaters,  as  implied,  are  coils  installed  within  the  supply  ventilation  ducting,  with 
separate  ventilation  fans  installed  upstream  of  coils.  Duct  heaters  provide  heat  to  the 
cutter’s  living  spaces. 

2.  Unit  Heaters 

a.  Dimensions 

Unit  heaters  are  manufactured  by  the  New  York  Blower  Company  of 
Willowbrook,  Illinois.  Two  models  are  employed:  Model  B-25  (quantity  two)  and  Model 
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SOLENOID  VALVE 


Figure  1.  Steam  Heating  System  Functional  Diagram  [Ref.  1] 


B-70  (quantity  five).  Dimensions  and  material  properties  of  xmit  heater  coils  and 
assemblies  are  shown  in  Appendix  B.  [Refs.  2-3] 
b.  Heating  Capacities 

Unit  heaters  are  sized  to  operate  with  steam  at  240  degrees  F  and  25  psig. 
For  each  unit  heater,  the  location,  air  flow  rate,  and  heat  transfer  rate  stated  on  the  ship’s 
prints  are  shown  in  Table  1 .  [Refs.  4-7] 


Heater 

Location 

Type 

Air  Flow 
(CFM) 

Heat  Transfer 
Stated  on  Prints 
(BTU/hr) 

2-8-0 

Paint  Locker 

B-70 

1,505 

5,294 

2-40-1 

Engine  Room 

B-70 

1,505 

64,860 

2-40-2 

Engine  Room 

B-70 

1,505 

64,860 

2-60-1 

Engine  Room 

B-70 

1,505 

64,860 

2-60-2 

Engine  Room 

B-70 

1,505 

64,860 

2-73-1 

Motor  Room 

B-25 

400 

33,120 

2-80-1 

B-25 

400 

7,747 

Table  1 .  Unit  Heater  Operating  Parameters 


3.  Duct  Heaters 

a.  Dimensions 

Duct  heaters  are  manufactured  by  two  manufacturers:  Colmac  Coil 
Manufacturing  Incorporated  of  Colville,  Washington  (quantity  four)  and  Carrier 
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Corporation  of  Farmington,  Connecticut  (quantity  one).  Duct  heater  dimensions  and 
material  properties  are  shown  in  Appendix  B.  [Refs.  8-9] 
b.  Heating  Capacities 

Duct  heaters  are  also  sized  to  operate  with  steam  at  240  degrees  F  and  25 
psig.  For  each  duct  heater,  the  supply  ventilation  system,  location,  air  flow  rate,  heat 
transfer  rate,  and  air  temperature  rise  stated  on  the  ship’s  prints  are  shown  in  Table  2. 
[Refs.  4-7] 


Heater 

System 

Location 

Air  Flow 
(CFM) 

Heat  Transfer 
Stated  on  Prints 
(BTU/hr) 

Air  Temp.  Rise 
Stated  on  Prints 
(°F) 

01-24-1 

S  01-24-1 

Fan  Room 

3,680 

51,600 

62  to  75 

01-50-0 

S  01-49-0 

Exhst.  Uptake 

750 

60,750 

-30  to  45 

01-25-1 

S  01-27-1 

Fan  Room 

1,400 

113,400 

-30  to  45 

2-25-1 

S  01-27-1 

Aux.  Mach.  1 

1,050 

45,360 

40  to  80 

2-16-1 

S  01-27-1 

Anchor  Gear 

350 

20,034 

45  to  98 

Table  2.  Duct  Heater  Operating  Parameters 


C.  HEAT  EXCHANGER  DESIGN  CONSIDERATIONS 

1.  Conventional  Heat  Exchangers 

The  design  of  a  heat  exchanger  involves  consideration  of  both  the  heat  transfer 
rates  between  the  fluids  and  the  mechanical  pumping  power  expended  to  overcome  fluid 
friction.  First  examining  heat  transfer  rates,  there  is  a  marked  difference  in  heat 
exchanger  performance  depending  on  whether  the  fluids  involved  are  gases  or  liquids. 
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The  convection  heat  transfer  coefficient,  h,  for  gases  is  generally  one  or  two  orders  of 
magnitude  less  than  that  for  liquids.  Correspondingly  a  gas  heat  exchanger’s  convective 
thermal  resistance,  R„  defined  as  the  inverse  of  the  product  of  h  and  heat  exchanger 
surface  area,  A^,  (i.e.  R,  =  l/hAj)  is  one  or  two  orders  of  magnitude  greater  than  that  for  a 
liquid  heat  exchanger.  It  is  therefore  evident  that  for  the  heat  transfer  rate  of  a  gas  heat 
exchanger  to  be  equivalent  to  that  of  a  liquid  heat  exchanger,  the  heat  transfer  surface 
area  for  a  gas  heat  exchanger  needs  to  be  much  larger  than  the  heat  transfer  surface  area 
for  a  liquid  heat  exchanger.  [Refs.  10-11] 

Examining  the  power  required  to  overcome  fluid  -  heat  exchanger  friction, 
conventional  (e.g.  concentric  tube  or  shell  and  tube)  heat  exchangers  operating  with  high 
density  fluids  have  lesser  frictional  losses  compared  to  heat  exchangers  operating  with 
low  density  fluids.  The  pumping  power  required  to  move  high  density  fluid  (e.g.  liquid) 
over  a  given  heat  exchanger  is  considerably  less  than  the  pumping  power  required  to 
move  low  density  fluid  over  the  same  heat  exchanger.  There  can  therefore  be  an 
equipment  operating  cost  benefit  to  using  a  liquid  heat  exchanger  system  rather  than  a  gas 
heat  exchanger  system. 

Heat  exchangers  where  at  least  one  fluid  is  required  to  be  a  gas  therefore  need  an 
improved  design  to  make  up  for  inherent  drawbacks.  The  heat  transfer  obtained  per  unit 
of  heat  exchanger  surface  area  can  be  increased  by  increasing  the  fluid  flow  velocity. 

This  however  is  not  a  desirable  method  to  increase  heat  transfer  since  the  friction  power 
expended  to  increase  fluid  flow  velocity  increases  by  as  much  as  the  cube  of  velocity. 
Minimizing  friction  power  leads  to  limiting  flow  velocities,  and  this  combined  with  the 
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relatively  low  thermal  conductivity  of  most  gases,  results  in  low  heat  transfer  rate  per  unit 
of  heat  exchanger  surface  area.  There  is  therefore,  in  a  conventional  heat  exchanger  using 
gases,  poor  heat  transfer  performance  at  low  flow  velocities  and  an  uneconomical, 
less  than  acceptable  trade-off  for  increasing  heat  transfer  per  unit  surface  area  by 
increasing  flow  velocities. 

2.  Compact  Heat  Exchangers 

The  development  of  compact  heat  exchangers  is  in  response  to  the  need  to  attain 
higher  heat  transfer  rates  with  minimum  space  and  power  requirements.  Large  but 
compact  surface  areas  are  a  typical  characteristic  of  gas  heat  exchangers.  Heat 
exchangers  used  on  the  WTGB  cutter  class  heating  system  are  one  style  of  compact  heat 
exchanger  that  incorporates  large  gas-side  surface  areas  with  dense  arrays  of  continuous 
fins. 

It  is  first  noted  that  compactness  itself  leads  to  high  performance.  A  compact 
surface  has  small  flow  passages  and  the  heat  transfer  coefficient,  h,  varies  as  a  negative 
power  of  the  flow  passage  size.  A  customary  expression  for  the  size  of  a  non-circular 
flow  passage  is  the  hydraulic  diameter,  equaling  four  times  the  cross-sectional  area 
divided  by  the  wetted  perimeter.  It  is  also  true  that  a  smaller  hydraulic  diameter  increases 
friction  power,  but  the  benefits  that  compacmess  has  on  the  heat  transfer  coefficient 
generally  outweigh  the  detrimental  influence  of  small  hydraulic  diameter  on  friction 
power.  [Ref.  10] 

In  addition  to  the  influence  of  small  hydraulic  diameter,  increases  in  heat 
exchanger  performance  can  be  obtained  by  any  modification  of  the  surface  geometry  that 
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results  in  a  higher  heat  transfer  coefficient  at  a  given  flow  velocity.  One  widely  accepted 
modification  is  use  of  extended  surfaces  or  fins  so  that  in  addition  to  providing  increased 
heat  transfer  surface  area,  the  interrupted  surface  prevents  thickening  boundary  layers 
from  reducing  heat  transfer.  Finned  surfaces  also  increase  friction  power  and  thermal 
resistance  due  to  conduction,  but  a  small  improvement  in  the  heat  transfer  coefficient  can 
more  than  offset  these  negative  factors.  [Ref.  10] 

Other  methods  of  obtaining  increased  performance  by  change  of  flow  surface 
geometry  include  the  use  of  curved,  corrugated,  or  wavy  passages,  in  which  boundary 
layer  separation  and  turbulence  (promoters  of  heat  transfer)  are  induced.  Such  surfaces 
are  incorporated  in  the  duct  heaters,  but  not  in  the  unit  heaters  being  studied. 

A  common  descriptor  of  compact  heat  exchangers  is  the  area  density,  a,  which  is 
the  ratio  of  heat  transfer  surface  area  to  heat  exchanger  volume.  A  conventional  cutoff 
for  labeling  a  heat  exchanger  as  compact  is  an  area  density  value  greater  than  700  square 
meters  per  cubic  meter  (or  213  square  feet  per  cubic  feet).  This  is  not  a  staunch  rule 
however,  as  many  heat  exchangers  have  been  grouped  into  the  compact  category  with 
lesser  area  densities  [Ref.  10]. 

Compact  heat  exchanger  designs  provide  the  benefits  of  high  heat  transfer  rate 
with  minimum  volume  and  thus  are  very  well  suited  for  duct  heater  applications.  In 
shipboard  systems  where  volume  savings  are  invariably  sought,  use  of  compact  heat 
exchangers  is  very  common. 
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III.  HEAT  EXCHANGER  ANALYSIS 


A.  STANDARD  PRESENTATION  OF  PERFORMANCE  DATA 
AND  THE  REYNOLDS  ANALOGY 

In  engineering  practice,  it  is  often  desirable  to  use  a  common  presentation  of 
performance  data  so  as  to  avoid  confusion  associated  with  many  arbitrarily  defined 
parameters.  In  the  study  of  heat  transfer  and  flow-friction,  commonality  in  data 
presentation  is  found  using  the  Reynolds  Analogy.  As  presented  by  Incropera  and  Dewitt 
[Ref  11],  and  also  Kays  and  Crawford  [Ref.  12],  with  certain  restrictions  (noted  shortly), 
relations  that  govern  velocity  boundary  layer  behavior  are  the  same  as  those  that  govern 
the  thermal  boundary  layer.  From  this  it  is  known  that  non-dimensional  friction  and  heat 
transfer  relations  for  a  particular  geometry  are  closely  related.  Specifically,  the  Reynolds 
number,  Re,  the  velocity  boundary  layer’s  friction  coefficient,  Cj,  and  the  heat  transfer 
boundary  layer’s  Nusselt  number,  Nu,  are  related  as  follows: 

where:  (2) 

F 

NH  =  -T 

k 


13 


Replacing  Nu  by  the  Stanton  number,  St  and  introducing  the  Prandtl  number,  Pr ; 


Nu, 

h 

St  =  - ^  (4) 

where:  St  = 

(5) 

ReiPr 

9Vc^ 

Pr  = 

P  ^ 

(6) 

k 

The  relation,  assuming  the  pressure  drop  in  the  flow  direction  is  zero  and  that  Pr=\  (true 
for  most  gases),  now  takes  the  form: 

-^  =  St  (7) 

2 

This  expression,  known  as  the  Reynolds  analogy,  relates  key  parameters  of  the 
velocity  and  thermal  boundary  layers.  The  accuracy  of  this  expression  depends  on  the 
noted  restrictions,  that  the  pressure  gradient  in  the  flow  direction  is  zero  and  the  Prandtl 
number  equals  one.  It  has  been  shown  that  the  analogy  may  be  applied  over  a  wide  range 
of  Prandtl  numbers  if  a  correction  is  added.  With  this  correction  arises  the  modified 
Reynolds,  or  Chilton-Colburn  analogy,  shown  as  follows: 

=  StPr^'^  =  0.6  <  Pr  <  60  (8) 

where  is  the  Colburn  j  factor.  For  laminar  flow,  the  modified  Reynolds  analogy  is 
again  only  appropriate  when  the  pressure  drop  in  the  flow  direction  is  zero.  In  tmbulent 
flow,  conditions  are  less  sensitive  to  the  effect  of  pressure  gradients  and  the  equation 
remains  valid  for  small  pressure  drops. 


14 


The  benefits  of  this  analogy  lie  in  the  ability  to  deduce  heat  transfer  information 
from  skin  friction  information  and  vice  versa.  A  wide  variety  of  compact  heat  exchanger 
performance  data  was  compiled  by  Kays  and  London  [Ref  10]  by  way  of  plots  of 
Colbumy  factor  and  friction  factor  versus  Reynolds  number.  The  scope  of  this  thesis 
includes  presentation  of  Colburn  j  factor  versus  Reynolds  number  and  does  not  include 
study  of  friction  factor.  In  using  this  standard  presentation,  future  correlation  of  friction 
factors  in  this  standardized  form  is  possible. 

B.  MANUFACTURER’S  DATA 


Plots  of  Colburn  j  factor  versus  Reynolds  number,  or  any  related  information; 
were  sought  from  the  each  of  the  heat  exchanger  manufacturers  with  no  success.  It  was 
evident  that  in  order  to  see  desired  heat  exchanger  performance  characteristics,  more 
commonly  available  information  such  as  that  shown  in  Figure  2  would  have  to  be  recast 
to  the  accepted  non-dimensional and  Re  forms. 


Air  Temperature  Rise  at  5  PSIG,  0°  EDB 

(WForWRAIum.  Rns) 

Face  Velocity,  SFPM 


Row 


Fin 

200 

400 

600 

800 

1000 

1200 

104 

49.7 

37.6 

31.8 

27.8 

25.3 

23.2 

106 

68.1 

51,1 

42.9 

37.7 

34.1 

31.5 

108 

85.2 

63.6 

52.6 

45.8 

41.2 

38.0 

110 

101.1 

74.9 

61.3 

53.4 

47.7 

43.1 

112 

114.9 

84.7 

69.7 

60.0 

53.1 

48.0 

206 

115.8 

90.8 

77.4 

69.0 

62.7 

58.1 

208 

138.5 

109,0 

93.1 

82.9 

74.9 

69.0 

210 

156.9 

124.9 

106.7 

94,3 

85.2 

78.6 

212 

171.5 

138.1 

118.1 

104.2 

93.9 

86.3 

Figure  2.  Typical  Manufacturer’s  Heat  Exchanger  Performance  Data 
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C.  ANALYSIS  MODEL 


The  heat  exchangers  studied  are  represented  in  Figure  3: 


where:  T,  =  Air  (cold  fluid)  inlet  temperature 

Tj  =  Air  (cold  fluid)  outlet  temperature 
T3  =  Steam  (hot  fluid)  inlet  temperature 
T4  =  Steam  (hot  fluid)  outlet  temperature 

D.  COMPUTATIONS  AND  RESULTS 


1.  Relation  Definitions 

Relations  used  in  the  steam  system  heat  exchanger  computations  are  defined  as 


follows: 


Mean  air  temperature: 


(T,  +  T^) 

Tm  =  -4-^  (9) 
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Air  mass  flow  rate: 


rh  .  =  p  .  *  V  . 

air  r  air  fr  air 

Air  side  heat  transfer  rate: 

Q  .  =  m  .  *c  * (r,  - ; 

^air  air  2 


Air  side  convection  coefficient: 

Qa 


h_._  = 


'air 


(Ts-T,)  *A^ 


Mass  velocity: 


G  = 


w. 


Reynolds  number: 


Re 


G  *D. 


Prandtl  number: 


Pr 


k 


Stanton  number: 


St 


h 

G  *c 

p 


Colburn  /  factor: 

jfj  =  St*Pr^'^ 


(10) 

\)  (11) 

(12) 

(13) 

(14) 

(15) 

(16) 

(17) 
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2.  Sequence  of  Computations 

Computations  for  the  steam  system  heat  exchanger  analyses  were  performed  with 
assistance  of  Fortran  computer  codes.  The  sequence  of  computations  in  the  Fortran  code 
titled  “hxair.f’,  shown  as  part  of  Appendix  C,  is  summarized  in  Figure  4.  Manufacturer’s 
data  for  each  heat  exchanger  similar  to  that  shown  in  Figure  2  were  input  to  the  program. 
Dimensional  characteristics  of  the  heat  exchangers  input  to  the  program  were  determined 
from  manufacturer’s  drawings  (see  Appendix  B  for  a  tabular  sxraimary  of  all  pertinent 
dimensions). 

3.  Tabular  and  Graphical  Results 

Program  output,  shown  in  tabular  form  in  Appendix  C,  included  the  Colburn  j 
factor  and  Reynolds  numbers  needed  for  further  analysis  of  a  hydronic  system  in  a 
subsequent  chapter.  Plots  of vs.  Re  for  each  heat  exchanger  studied  are  shown  in 
Figures  5  through  11. 

4.  Comparison  with  Established  Results 

Results  obtained  from  the  steam  system  analysis  were  compared  with  the 
established  results  of  Kays  and  London  [Ref  10].  For  a  similarly  configured  heat 
exchanger,  a  Colburn j  factor  versus  Reynolds  number  plot  is  shoivn  in  Figure  12.  It  is 
evident  that  orders  of  magnitude  of  Colburn  J  factors  at  respective  Reynolds  numbers  and 
overall  trends  of  data  for  this  analysis  compare  favorably  with  previously  established 
results. 
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Figure  4.  Steam  System  Computations  Flow  Chart 
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COLBURN  j  FACTOR  vs.  REYNOLDS  NUMBER 


Figure  6.  Colburn  j  Factor  vs.  Reynolds  Number  -  Unit  Heater  B70 


COLBURN  j  FACTOR  vs.  REYNOLDS  NUMBER 


Figure  7.  Colburn  j  Factor  vs.  Reynolds  Number  -  Duct  Heater  01-24-1 


COLBURN  j  FACTOR  vs.  REYNOLDS  NUMBER 


Figure  8.  Colburn  j  Factor  vs.  Reynolds  Number  -  Duct  Heater  01-25-1 


COLBURN  j  FACTOR  vs.  REYNOLDS  NUMBER 


Figure  9.  Colburn  j  Factor  vs.  Reynolds  Number  -  Duct  Heater  01-50-0 


COLBURN  j  FACTOR  vs.  REYNOLDS  NUMBER 


Figure  10.  Colburn  j  Factor  vs.  Reynolds  Number  -  Duct  Heater  2-25-1 


Finned  circular  tubes,  surface  8.0-3/8T. 


Tube  outside  diameter  =  0.402  in  =  1 0.2  x  10'^  m 
Fin  pitch  =  8.0  per  in  =  31 5  per  m 

Flow  passage  hydraulic  diameter,  =  0.01 192  ft  =  3.632  x  10"^ m 
Fin  thickness  =  0.01 3  in  =  0.33  x  1 0"^ m 
Free  flow  area/frontal  area,  a  =  0.534 

Heat  transfer  area/total  volume,  a  =  179  ft^/ft^  =  587  m^/m^ 

Fin  area/total  area  =  0.913 

Note;  Minimum  free-flow  area  in  spaces  transverse  to  flow. 


Figure  12.  Established  Colburn  j  Factor  vs.  Reynolds  Number  Relation  [After  Ref.  10] 
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IV.  HYDRONIC  SYSTEM  ANALYSIS 


A.  APPROACH 

To  analyze  the  compartment  heating  system  as  a  hydronic  (circulating  water) 
system,  the  Effectiveness  -  NTU  Method  is  utilized.  Required  for  this  method  are 
determinations  of  heat  exchanger  overall  heat  transfer  coefficients,  which  entails 
computing  other  parameters  such  as  air  and  water  convection  coefficients,  fm 
efficiencies,  wall  conduction  resistances,  and  fouling  resistances.  Results  from  the 
previous  chapter  are  an  integral  part  of  this  analysis,  and  are  used  specifically  to 
determine  air  side  convection  coefficients  at  specified  fan  operating  conditions. 

B.  OVERALL  HEAT  TRANSFER  COEFFICIENT 


1.  Newton’s  Law  of  Cooling 

A  standard  measure  of  a  heat  exchanger’s  total  thermal  resistance  to  heat  transfer 
between  two  fluids  separated  by  a  heat  exchanger  structure  is  the  overall  heat  transfer 
coefficient,  U.  In  an  expression  analogous  to  Newton’s  Law  of  Cooling,  the  heat  transfer 
between  two  fluids  separated  by  one  or  more  thermal  resistances,  R„  is: 


Q  - 


At 


=  UA  At 

s 


(18) 


The  result  apparent  from  the  above  expression,  that  is  applicable  to 

clean,  unfinned  surfaces  only,  but  can  be  used  as  the  basis  for  determining  a  heat 
exchanger’s  overall  heat  transfer  coefficient  with  additional  factors  considered  such  as  the 
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effects  of  fins  on  the  air  side  surfaces,  the  presence  of  fouling  on  both  air  and  water  sides, 
and  tube  wall  conduction. 

2.  Overall  Coefficient  Elements 

a.  Air  Side  Convection  Coefficients 

The  air  side  convection  coefficients,  are  determined  using  the 
Reynolds  number  vs.  Colburn  j  Factor  plots,  generated  as  part  of  the  steam  system 
analysis  in  the  previous  chapter.  The  utility  of  presenting  data  in  non-dimensional  form  is 
now  very  apparent,  as  the  performance  characteristics  of  each  heat  exchanger  formulated 
with  steam  now  is  relevant  to  the  hydronic  analysis. 

An  actual  operating  condition  of  the  fan  supplying  air  to  a  particular  heat 
exchanger,  in  the  form  of  a  quantity  of  cubic  feet  per  minute  of  air  at  a  specified  air 
temperature  (taken  from  manufacture’s  data)  is  used  to  compute  the  air  flow’s  mass  flow 
rate,  rh  .  ,  then  mass  velocity,  G,  then  Reynolds  number.  Re.  Entering  the  vs.  graph 


with  the  computed  Reynolds  number,  a  value  for is  obtained,  fi:om  which  is 
computed  using  the  Stanton  number  and  Colburn  j  factor  relations,  simplified  into  the 
following; 


h 


air 


Prin 


(19) 


Computations  of  the  air  side  convection  coefficient  for  each  of  the 
system’s  heat  exchangers  are  shown  in  Table  3.  Note  that  model  B-25  unit  heaters  (i.e.  2- 
8-0  and  2-73-1)  and  model  B-70  unit  heaters  (2-80-1, 2-40-1, 2-40-2, 2-60-1, 2-60-2) 
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each  have  one  type  fan  and  thus  a  common  volumetric  flow  rate  of  air  (“CFM  of  Air  @ 
70  °F”),  and  were  each  grouped  into  one  row  in  the  table. 


Heater 

Number 

CFM  of  air 
@70°F 

rh 

(Ibm/hr) 

G 

(Ibm/ft^-sec) 

Re 

7h 

Kir 

(BTU/hr-tf-R) 

Type  B-25* 

400 

1,798 

1.16 

2,981 

0.0048 

6.0 

Type  B-70** 

1,505 

6,763 

2.25 

5,795 

0.0037 

9.0 

01-24-1 

3,680 

16,538 

1.63 

2,762 

0.0040 

7.0 

01-25-1 

1,400 

6,292 

1.62 

1,691 

0.0047 

8.2 

01-50-0 

750 

3,371 

1.16 

1,211 

0.0056 

7.0 

2-25-1 

1,050 

4,719 

1.73 

3,654 

0.0054 

10.1 

2-16-1 

350 

1,573 

0.82 

851 

0.0062 

5.4 

Type  B-25  includes  unit  heaters  2-8-0  and  2-73-1 . 

Type  B-70  includes  unit  heaters  2-80-1, 2-40-1, 2-40-2,  2-60-1, 2-60-2. 

Table  3.  Air  Side  Convection  Coefficient  Computations 


b.  Water  Side  Convection  Coefficients 
The  water  side  convection  coefficients,  h„„.,  are  computed  based  on 
internal  flow  relations  presented  by  Incropera  and  DeWitt  [Ref.  11].  Initial  computations 
of  water  side  convection  coefficients  are  based  on  a  manufacturer’s  recommended 
maximum  water  flow  rate  of  ten  gallon  per  minute.  These  initial  water  side  convection 
coefficients  will  be  modified  in  some  cases  based  on  an  optimized  water  mass  flow  rate, 
as  will  be  shown  in  a  later  section.  It  is  also  noted,  as  it  was  in  the  previous  chapter,  that 
the  water  side  convection  coefficient  is  generally  one  or  two  orders  of  magnitude  higher 
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than  the  gas  side  convection  coefficient.  The  water  side  convection  thermal  resistance  (R, 
=  1/hA)  will  therefore  be  one  or  two  orders  of  magnitude  lower  than  that  of  the  air  side, 
and  a  relatively  small  contributor  to  the  sum  of  thermal  resistances  used  to  compute  the 
overall  heat  transfer  coefficient,  U,  as  will  be  shown  in  a  following  section. 

(1)  Reynolds  Number  for  Internal  Flow.  The  Reynolds  number  for 
internal  flows  can  be  computed  using  the  following  relation,  where  the  hydraulic 
diameter,  D/,  —  4A/P. 


4m 


wtr 


%D,ix 


h^wtr 


(20) 


(2)  Nusselt  Number.  The  Nusselt  Number,  Nup,  relations  for 
internal  laminar  flow  <  2,300)  and  internal  turbulent  flow 

(Rcj^  >  2,300  -  Gnielinski  Correlation)  are  then  Utilized  to  solve  for 


Nuj^  =  (21) 

k 

=  4.36 
=  6.49 

(/78)(Reo  -  1000)Pr 
1  +  12.7(/78)^'^(/’r^^^  -  1) 

where:  /=  (0.791n/?e^  -  1.64)'^ 


Re^  <  2300,  circular  duct  heater  tubes 
Reo  <  2300,  rectangular  unit  heater  tubes 

Re^  >  2300,  duct  heater  and  unit  heater  tub^ 
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(3)  Sequence  of  Computations.  Computations  of  the  water  side 
convection  coefficients  for  each  of  the  system’s  heat  exchangers  were  completed  with 
assistance  of  Fortran  computer  codes.  The  sequence  of  computations  in  the  Fortran  code 
titled  “hxwtr.f’j  shown  as  part  of  Appendix  C,  is  summarized  in  Figure  13.  Dimensional 
characteristics  of  the  heat  exchangers  input  to  the  program  were  determined  from 
manufacturer’s  drawings  (see  Appendix  B  for  a  tabular  summary  of  all  pertinent 
dimensions). 

(4)  Tabular  Results.  Excerpts  from  program  output  is  shown  in 
Tables  4  through  6.  Note  that  the  water  tubing  dimensions  for  type  B-25  and  B-70  unit 
heaters  are  identical,  and  thus  these  results  are  grouped  into  Table  4.  Similarly,  tubing 
dimensions  for  all  duct  heaters  with  the  exception  of  heater  01-24-1  are  identical  and 
these  results  are  grouped  into  Table  5.  Results  for  heater  01-24-1,  with  its  own  unique 
tubing  dimensions  are  shown  in  Table  6.  Complete  program  output  is  shown  in 
Appendix  C. 

c.  Extended  Surfaces  (Fins) 

As  discussed  in  the  previous  chapter,  extended  surfaces  or  fins  are  often 
added  to  heat  exchanger  surfaces  to  decrease  thermal  resistance  to  convection  heat  trans¬ 
fer.  A  representative  fin  array  for  the  heat  exchangers  being  analyzed  is  shown  in  Figure 
14.  Expressions  for  fin  efficiencies  are  as  presented  by  Incropera  and  DeWitt  [Ref  11]. 
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Flow  Rates 


Mass 

(Ibm/hr) 

Volume 

(gal/min) 

Rcd 

Nud 

(BTU/hr-tf-R) 

48.8 

0.100 

1482.9 

6.49 

53.4 

97.5 

0.200 

2965.8 

15.13 

124.4 

146.3 

0.300 

4448.7 

23.88 

196.4 

780.2 

1.600 

23726.5 

107.45 

883.6 

828.9 

1.700 

25209.4 

113.04 

929.6 

4876.0 

10.00 

148290.7 

492.92 

4053.2 

Table  4.  Unit  Heater  Water  Side  Convection  Coefficient  Computations 


Flow  Rates 

Mass  Volume 

flbni/hrl  teal/minl 

Ren 

Nud 

(BTU/hr-tf-R) 

48.8 

0.100 

1505.1 

4.36 

36.4 

97.5 

0.200 

3010.2 

15.41 

128.6 

146.3 

0.300 

4515.3 

24.25 

202.4 

780.2 

1.600 

24081.5 

108.80 

908.0 

828.9 

1.700 

25586.6 

114.46 

955.3 

4876.0 

10.00 

150509.5 

499.06 

4165.2 

Table  5.  Duct  Heater  Water  Side  Convection  Coefficient  Computations  (with  exception 

of  heater  01-24-1) 
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Flow  Rates 

Mass  Volume 

(Ibm/hr)  (gal/min) 

Rcj) 

Nud 

hwtr 

(BTU/hr-ft^-R) 

48.8 

0.100 

807.1 

4.36 

19.5 

97.5 

0.200 

1614.2 

4.36 

19.5 

146.3 

0.300 

2421.2 

11.55 

51.7 

780.2 

1.600 

12913.3 

64.15 

287.1 

828.9 

1.700 

13720.4 

67.58 

302.5 

4876.0 

10.00 

80708.0 

297.17 

1329.9 

Table  6.  Duct  Heater  01-24-1  Water  Side  Convection  Coefficient  Computations 


(1)  Single  Fin  Efficiency.  The  single  fin  efficiency,  7^,  is  given  by  the  expression: 


tanh  mL 


n.  = 


mL 


'f"‘c 


(22) 


where: 


m  = 


C23) 


+ 


fix 


fix 


(24) 
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Ac  =  wt 


Figure  14.  Representative  Fin  Array 


(2)  Fin  Overall  Surface  Efficiency.  Characterizing  the 
performance  of  an  array  of  fins  is  done  with  the  expression  fm  overall  surface  efficiency, 
where; 


^max 


hA  (T.  -  T) 

s  .  ^  base 


(25) 


Breaking  down  the  total  heat  transfer  into  contributions  from  the 
finned  and  unfinned  (base)  surfaces  and  noting  that  =  A^^  '^■^base >  expression  can 


be  further  simplified  yielding: 


n. 


1  - 


(1  -^f) 


(26) 


(3)  Sequence  of  Computations.  Computations  of  single  and 
overall  fm  efficiencies  for  each  of  the  system’s  heat  exchangers  were  completed  with 
assistance  of  Fortran  computer  codes.  The  sequence  of  computations  in  the  Fortran  code 
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titled  “hxU.f’,  shown  as  part  of  Appendix  C,  is  summarized  in  Figure  15.  Dimensional 
characteristics  of  the  heat  exchangers  input  to  the  program  were  determined  from 
manufacturer’s  drawings  (see  Appendix  B  for  a  tabular  summary  of  all  pertinent 
dimensions). 

(4)  Tabular  Results.  A  summary  of  program  output  for  each  of  the 
system’s  heat  exchangers  is  shown  in  Table  7.  Efficiencies  were  computed  using 
operating  conditions  and  air  side  convection  coefficients  of  Table  3.  Complete  program 
output  is  shown  in  Appendix  C. 


Heater 

B-25*  Unit  Heaters 

0.87 

0.89 

B-70**  Unit  Heaters 

0.82 

0.84 

01-24-1  Duct  Heater 

0.92 

0.93 

01-25-1  Duct  Heater 

0.90 

0.90 

01-50-0  Duct  Heater 

0.91 

0.92 

2-25-1  Duct  Heater 

0.92 

0.92 

2-16-1  Duct  Heater 

0.93 

0.93 

Table  7.  Heater  Single  Fin  and  Overall  Fin  Efficiencies 


d.  Fouling  Resistance 

During  normal  heat  exchanger  operation,  surfaces  are  often  subject  to 
fouling  by  fluid  impurities,  rust  or  scale  formation,  and  other  reactions  between  the  fluid 
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and  the  wall  material.  Depositions  sueh  as  these  can  greatly  increase  the  resistance  to 
heat  transfer  between  the  fluids.  This  effect  can  be  treated  by  introducing  an  additional 
thermal  resistance,  termed  the  fouling  factor,  Ry,  the  value  of  which  depends  on  the 
operating  temperature,  fluid  velocity,  and  length  of  service  of  the  heat  exchanger. 

Fouling  factors  recommended  by  the  Tubular  Exchanger  Manufacturer’s  Association  are 
shown  below  in  Table  8.  The  shaded  fouling  factors  are  those  used  for  the  heating  system 
analyzed  in  this  thesis. 

e.  Tube  Wall  Conduction 

Derivations  of  the  tube  wall  conduction  resistance,  by  Incropera  and 


Fluid 

Fouling  Factor  (hr-ft^-°F/BTU) 

Below  125  °F  Above  125  °F 

Seawater 

0.0005 

0.001 

DistUkdW^r 

■■■■■■ 

Treated  Boiler  Feedwater 

0.001 

0.001 

City  or  Well  Water 

0.001 

0.002 

Hard  Water 

0.003 

0.005 

Ak  ' 

0.002 

p  M2  [ 

Diesel  Exhaust 

0.01 

0.01 

Clean  Steam 

0.0005 

0.0005 

Table  8.  Typical  Fouling  Factors 
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DeWitt  [Ref.  11]  produce  the  following  expressions  for  duct  heater  tubes  (cylindrical) 
and  unit  heater  tubes  (modeled  as  plane  wall). 


R 


wall 


(27) 


(cylindrical  tube  wall  conduction  resistance) 


=  ( 77)  (2«) 

\  KA^/  tube 

(plane  wall  tube  conduction  resistance) 

3.  Summation  of  Elements 

The  summation  of  thermal  resistances  including  effects  of  air  side  convection 
with  fins,  water  side  convection  without  fins,  fouling  on  both  air  and  water  sides,  and 
tube  wall  conduction  yields  the  following: 


1 


UA 


U_,A 


U  ,  A 

wtr 


R. 


(T]  hA  )  .  (r\  A  ) 

^  *0  S' air  ^  *0  s' air 


R^  1 

+  R  ^  +  - - 

wall  A 


(29) 


The  calculation  of  the  UA^  product  does  not  require  designation  of  the  air  side  or 
water  side,  i.e.  UA  -  U  A  -  U  ,  A  .  However,  the  calculation  of  an  overall 

^  ^air  ^wtr 


coefficient,  U,  does  depend  on  whether  it  is  based  on  the  air  side  or  the  water  side  surface 
area  since  ^  ^  j  ^  .As  mentioned  previously,  the  convection 

air  wtr  J  s  .  s  . 


air  wtr 
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coefficient  for  the  air  side  is  generally  much  smaller  than  that  for  the  water  side,  and  thus 
dominates  determination  of  the  overall  coefficient. 

4.  Overall  Coefficient  Computations  and  Results 

a.  Sequence  of  Computations 

Computations  of  the  overall  coefficients  for  each  of  the  system’s  heat 
exchangers  were  completed  with  assistance  of  Fortran  computer  codes.  The  sequence  of 
computations  in  the  Fortran  code  titled  “hxU.f shown  as  part  of  Appendix  C,  is  summa¬ 
rized  in  Figure  15.  Program  inputs  are  the  previously  computed  values  of  h^;,.,  tabu¬ 
lated  values  of/? ,  and/?,  ,  and  dimensional  and  material  characteristics  (see  Appendix 

J  air  J  wtr 

B  for  a  tabular  summary  of  all  pertinent  dimensions). 

b.  Tabular  Results 

A  summary  of  program  output  for  each  of  the  system’s  heat  exchangers  is 
shown  in  Table  9.  Complete  program  output  is  shown  in  Appendix  C. 

C.  THE  EFFECTIVENESS  -  NTU  METHOD 

1.  Applicability 

The  Effectiveness  -  NTU  analysis  method  is  used  to  take  the  overall  heat  transfer 
coefficient  parameter  and  determine  heat  exchanger  heat  transfer  performance 
characteristics.  This  method  is  chosen  over  the  log-mean  temperature  difference  method 
based  on  a  detailed  comparison  between  the  two  methods  presented  by  Kays  and  London 
[Ref.  10].  In  the  comparison,  steps  for  each  method  are  evaluated  for  a  variety  of 
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Figure  15.  Fin  Efficiency  and  Overall  Coefficient  Computations  Flow  Chart 
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Heater 

(BTU/hr-ft^-R) 

B-25  Unit  Heaters* 

4.9 

B-70  Unit  Heaters** 

6.7 

01-24-1  Duct  Heater 

5.7 

01-25-1  Duct  Heater 

6.4 

01-50-0  Duct  Heater 

5.6 

2-25-1  Duct  Heater 

8.3 

2-16-1  Duct  Heater 

4.5 

Type  B-25  include  unit  heaters  2-8-0  and  2-73-1 . 

Type  B-70  include  unit  heaters  2-80-1, 2-40-1, 2-40-2, 2-60-1, 2-60-2. 


Table  9.  Heat  Exchanger  Overall  Coefficient  Results 


analysis  starting  points.  For  this  analysis,  where  only  the  inlet  hot  and  cold  fluid 
temperatures  are  known,  the  Effectiveness  -  NTU  method  is  more  straightforward, 
eliminating  a  tedious,  iterative  procedure  where  outlet  temperatures  would  be  estimated 
and  then  adjusted  until  the  heat  transfer  rate  corresponds  to  the  inlet/outlet  temperature 
difference.  In  the  Effectiveness  -  NTU  (abbreviation  for  “number  of  transfer  units”) 
approach,  from  knowledge  of  the  heat  exchanger  type,  size,  and  fluid  flow  rate,  the  NTU, 
maximum  heat  transfer  rate,  n  »  ^nd  effectiveness,  e  (both  defined  shortly)  are  used  to 

^max 

determine  the  actual  heat  transfer  rate. 
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2. 


Maximum  Possible  Heat  Transfer  Rate 


The  maximum  possible  heat  transfer  rate  for  a  given  heat  exchanger  is  achieved 
when,  by  definition,  the  entering  cold  air  temperature  is  raised  to  equal  its  highest 
possible  value,  i.e.,  that  of  the  entering  hot  water  temperature.  Referring  to  Figure  3, 
entering  air  would  rise  in  temperature  by  the  quantity  (T^  -  .  Using  a  relationship  for 

the  heat  transfer  rate  similar  to  that  presented  earlier: 

Q  =  mc^(T^  -  r,)  (30) 

where  the  term  mc^  is  defined  as  the  heat  capacity  rate,  C,  and  will  be  either  based  on 
the  mass  flow  rate,  m ,  and  specific  heat,  ,  for  air  or  for  water.  As  developed  by 
Incropera  and  DeWitt  [Ref.  11],  the  maximum  heat  transfer  rate,  q  ,  is  then: 

^max 

Q  ^  C  .  (T-  T.)  (31) 

^max  min'"  3  I'  ^  ^ 

where  is  equal  to  either  or  whichever  is  smaller. 

3.  Effectiveness 

The  effectiveness,  e,  is  defined  as  the  ratio  of  the  actual  heat  transfer  rate  for  a 
heat  exchanger  to  the  maximum  possible  heat  transfer  rate. 

e  =  (32) 

Q 

^max 

The  effectiveness,  which  is  dimensionless,  must  be  in  the  range  0  <  e  <  1 .  It  is 
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useful  because,  if  e,  T^,  and  are  known,  the  actual  heat  transfer  rate  may  be 

determined  from  the  expression: 

Q  =  eO  =  eC  .  (T^  -  T,)  (33) 

4.  Number  of  Transfer  Units 

The  number  of  transfer  units,  NTU,  is  a  dimensionless  parameter  that  is  widely 
used  for  heat  exchanger  analysis  and  is  defined  as: 

UA 

NTU  =  - ^  (34) 

C  . 

min 


Effectiveness  -  NTU  Relations 


The  Effectiveness  -  NTU  relationship  for  the  type  heat  exchangers  being  analyzed, 
i.e.  a  single  pass  cross  flow  heat  exchanger  with  both  fluids  unmixed  is: 


e  =  1  -  exp 


(TVTC/f^^jexp  -C^iNTUf^  -  l} 


(35) 


where: 


(36) 


The  NTU,  e,  relation  may  be  represented  graphically  as  shown  in  Figure  16. 

6.  Computations  and  Results 

a.  Sequence  of  Computations 

Computations  involving  the  Effectiveness  -  NTU  analysis  of  the  system’s 


45 


heat  exchangers  were  completed  utilizing  a  spreadsheet  format,  shown  in  Appendix  D. 
Inputs  are  heat  exchanger  dimensions,  water  and  air  inlet  temperatures  T,  and  Tj 
respectively,  air  and  water  volumetric  flow  rates,  and  previously  determined  values  of  jf,. 


Figure  16.  NTU,  e,  Q.  Relation  of  a  Single-Pass  Cross  Flow  Heat  Exchanger  with  Both 

Fluids  Unmixed 

and  Using  relations  defined  in  Sections  IV-C-2  through  IV-C-5,  NTU- 
effectiveness  parameters  were  computed  and  finally  the  heat  transfer  rate  and  outlet  air 
and  water  temperatures  Tj  and  T4  were  determined. 
b.  Tabular  Results 

Effectiveness  -  NTU  analysis  results  for  the  hydronic  system  are  shown  in 

Table  10. 
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V.  SYSTEM  COMPAMSONS  AND  DISCUSSION 

A.  OUTLINE 

The  hydronic  system  heat  exchanger  performance  data  generated  in  the  previous 
chapter  will  be  compared  with  steam  system  performance  parameters.  Specifically,  the 
heat  rate  outputs  (Q)  and  the  final  air  temperatures  (T^)  of  the  hydronic  system  will  be 

compared  with  those  of  the  steam  system.  In  cases  where  the  hydronic  system 
performance  matches  or  exceeds  steam  system  performance  parameters  stated  on  the 
ship’s  prints,  water  flow  rates  will  be  optimized.  In  instances  where  the  hydronic  system 
performance  falls  short  of  that  of  the  steam  system,  further  analysis  will  be  conducted  to 
attempt  to  increase  hydronic  system  heat  transfer. 

B.  FURTHER  EXAMINATION  OF  STEAM  SYSTEM 
CAPABILITIES 

1.  Modifications  to  Effectiveness-NTU  Method 

Steam  system  heat  exchanger  heat  transfer  performance  data  is  essential  in  order 
to  assess  the  capabilities  of  the  hydronic  system.  To  compute  steam  system  performance 
parameters,  the  Effectiveness-NTU  analysis  method  is  again  used,  modified  for  steam  as 
follows: 

1 .  The  convection  coefficient  for  steam  is  taken  to  be  very  large  (on  the  order  of 
10'*),  thereby  decreasing  the  air  side  thermal  resistance,  contributing  negligibly 
to  the  overall  heat  transfer  coefficient. 
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2.  The  specific  heat  for  a  substance  undergoing  a  phase  change,  such  as  steam 
condensing,  is  infinite.  The  heat  capacity  rate  for  steam,  is  therefore 
also  infinite,  yielding  a  heat  capacity  ratio,  C^,  effectively  equal  to  zero. 

3.  In  the  case  of  ®  ,  the  effectiveness,  e  =  1  -  exp  ( -NTU ) . 

4.  The  inlet  air  and  steam  temperatures  are  as  stated  previously  in  Chapter  II. 

5.  The  steam  mass  flow  rate,  ^  ^  =  qiQi^  ) ,  where  hf^  is  the  enthalpy 

difference  due  to  condensation  from  saturated  steam  to  saturated  water. 

2.  Sequence  of  Computations 

Computations  involving  the  Effectiveness-NTU  analysis  of  the  system’s  heat 
exchangers  were  completed  utilizing  a  spreadsheet  format  shown  in  Appendix  D.  Inputs 
are  heat  exchanger  dimensions,  steam  and  air  inlet  temperatures,  Tj  and  Tj  respectively, 
as  well  as  air  volumetric  flow  rates,  and  previously  determined  values  of  and 

Using  relations  defined  in  Chapter  III,  Effectiveness-NTU  parameters  were  computed 
and  finally  the  heat  transfer  rate  and  outlet  air  temperature  for  each  heat  exchanger 
were  determined. 

3.  Tabular  Results 

Effectiveness-NTU  analysis  results  for  the  steam  system  are  shown  in 
Table  11. 
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C.  COMPARISONS  OF  HYDRONIC  AND  STEAM  SYSTEMS 


1.  Tabular  Results 

Comparisons  of  the  hydronic  and  steam  systems’  heat  transfer  performances  are 
shown  in  Tables  12  and  13. 


Heater 

^hydronic 

(BTU/hr) 

6 

^  Steam 

(BTU/hr) 

Q  hydronic 

Q 

^  steam 

B-25 

11,713 

16,722 

0.70 

B-70 

28,214 

40,445 

0.70 

Table  12:  Unit  Heater  Hydronic  vs.  Steam  Performance  Data 


Heater 

^hydronic 

(BTU/hr) 

Ar 

^^^hydronic 

(°F) 

6 

^  steam 

(BTU/hr) 

At 

(°F) 

^hydronic 

^  steam 

At 

hydronic 

At  . 

“"■sleam 

01-24-1 

74,778 

19 

121,615 

31 

0.61 

0.61 

01-25-1 

66,393 

44 

87,498 

58 

0.76 

0.76 

01-50-0 

43,215 

53 

56,155 

69 

0.77 

0.77 

2-25-1 

22,526 

20 

31,449 

28 

0.71 

0.71 

2-16-1 

15,152 

40 

21,428 

57 

0.70 

0.70 

Table  13:  Duct  Heater  Hydronic  vs.  Steam  Performance  Data 
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2.  Comparison  Trends 

Comparisons  between  hydronic  and  steam  system  results  bear  some  consistencies. 
Notably,  the  ratio  of  the  hydronic  to  steam  heat  performance  values  average 
approximately  70  percent.  This  is  a  performance  difference  which  warrants  examination 
to  determine  what  the  possibilities  are  for  improved  performances  and  what  are  the  causes 
for  the  shortfalls. 

D.  IMPROVING  PERFORMANCE  OF  HYDRONIC  SYSTEM 

1.  Altering  Air  Mass  Flow  Rate 

a.  Analysis  Method 

The  Effectiveness-NTU  analysis  makes  very  apparent  that  the  air-water 
heat  exchangers  studied  are  “air  side  limited”.  This  can  be  explained  qualitatively  by 
examining,  for  given  heat  exchanger  dimensions  and  flow  conditions,  the  relative 
magnitudes  of  the  contributors  to  the  basic  relations  used  in  the  analysis.  The  inverse 
relationships  of  thermal  resistances  in  the  computation  of  overall  coefficient,  U,  heavily 
skews  the  resulting  value  of  U  to  the  lower  convection  coefficient  in  the  equation ,  in  this 
case  The  number  of  transfer  units,  NTU,  is  directly  proportional  to  U,  thus  a  lower 

value  of  U  results  in  a  lower  NTU.  Examining  the  graphical  representation  of  the  NTU-e- 
Q  relation  (see  Figure  16),  lower  values  of  NTU  yield  lower  values  of  effectiveness  e. 
The  heat  transfer  rate,  q  ,  is  then  the  product  of  n  and  e. 

^  ^max 

A  means  to  increase  the  air  side  convection  coefficient,  and  thereby 
increasing  the  overall  coefficient  and  the  heat  transfer  rate,  is  to  increase  the  rate  of  air 


53 


flow  passing  over  the  heat  exchanger.  A  change  in  the  air  side  convection  coefficient 
with  increasing  air  flow  rate  can  be  predicted  by  examining  the  plots  of  Colburn  j  factor 
vs.  Reynolds  number  and  the  relation  from  which  is  derived,  i.e. 

^air  ~  negative  slope  of  the  plots  indicate  that  jfj  decreases  with 

increasing  Reynolds  number,  but  it  is  also  true  that  the  mass  velocity,  G,  increases  with 
increasing  Reynolds  number.  The  greatest  increases  in  will  be  realized  when  Re  and 
therefore  G  increase  with  minimal  decrease  in  jff.  This  leads  to  the  qualitative  deduction 
that  increases  in  are  more  profound  if  the  plot  of  Colburn  J  factor  vs.  Reynolds 
number  is  nearer  to  horizontal. 

b.  Analysis  Results 

The  net  effects  on  of  the  coxmteracting  trends  of  and  G,  and  the 

effects  of  the  plot’s  slope  can  best  be  seen  with  sample  computations.  Results  of 
Effectiveness-NTU  analysis  computations  at  air  flow  rates  of  up  to  twice  the  original  flow 
rates  for  duct  heater  B-25  are  shown  in  Table  14.  Computations  for  the  remaining  heat 
exchangers  are  shown  in  tables  in  Appendix  D.  Plots  of  heat  transfer  rate  vs.  air  flow  rate 
for  each  heat  exchanger  are  shown  in  Figures  1 7  through  23. 

c.  Air  Temperature  Rise  Considerations 

The  increased  heat  transfer  rate  at  increased  air  flow  rate  produces  a 
decrease  in  the  rise  in  air  temperature.  This  can  be  explained  by  examining  the 
relationship  between  the  heat  transfer  rate  and  the  change  in  air  temperature  presented  in 
a  previous  chapter,  i.e  t^=T^+  {QIC  •  Th®  transfer  rate  ,  q  ,  increases  but  the 
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Heat  Transfer  Rate  (BTU/hr) 


HEAT  TRANSFER  RATE  vs.  AIR  FLOW  RATE 
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Figure  17.  Heat  Transfer  Rate  vs.  Air  Flow  Rate  -  Unit  Heater  B25 
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Figure  18.  Heat  Transfer  Rate  vs.  Air  Flow  Rate  -  Unit  Heater  B70 
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Figure  19.  Heat  Transfer  Rate  vs.  Air  Flow  Rate  -  Duct  Heater  01-24-1 
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Figure  20.  Heat  Transfer  Rate  vs.  Air  Flow  Rate  -  Duct  Heater  01-25-1 
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HEAT  TRANSFER  RATE  vs.  AIR  FLOW  RATE 
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Figure  22.  Heat  Transfer  Rate  vs.  Air  Flow  Rate  -  Duct  Heater  2-25-1 
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Figure  23.  Heat  Transfer  Rate  vs.  Air  Flow  Rate  -  Duct  Heater  2-16-1 
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air’s  heat  capacity,  C^.^,  also  increases  with  larger  air  flow  rates,  resulting  in  a  net 
reduction  in  air  temperature  rise. 

The  reason  for  providing  heat  to  a  particular  volume  of  air  needs  to  be 
considered  to  determine  whether  a  decrease  in  is  acceptable  at  a  higher  q.  If  the 

purpose  of  providing  heat  to  a  volume  of  air  is  solely  to  increase  the  air  temperature  of 
potentially  very  cold  air,  as  is  the  case  with  duct  preheaters,  then  the  decrease  in  at  a 

higher  q  would  not  be  acceptable.  If  the  purpose  of  providing  heat  to  a  space  is  to 

provide  incident  heat  to  surfaces  in  the  space,  with  the  air  temperature  already  at 
reasonably  comfortable  level,  (e.g.  60 °F),  as  could  be  the  case  with  unit  heaters,  then  the 
decrease  in  Ar^.^  at  a  higher  q  would  be  acceptable.  In  the  latter  case,  the  air  provided  to 

the  space  at  an  increased  q  .has  increased  enthalpy,  thereby  able  to  increase  the  enthalpy 

of  a  volume  of  ‘unheated’  air  to  a  proportionally  higher  value. 
d.  Practicality 

Increasing  the  flow  rate  of  air  passing  over  the  duct  heaters,  if  involving 
the  replacement  of  existing  fans  with  larger  fans,  is  not  seemingly  a  practical  endeavor. 
The  expense,  effort,  and  potential  ramifications  of  duct  fan  replacement  make  this  course 
of  action,  at  an  initial  evaluation,  not  viable. 

Increasing  the  flow  rate  of  air  passing  over  the  unit  heaters,  due  to  much 
better  accessibility  and  smaller  scale  of  the  project,  is  seemingly  a  practical  endeavor.  An 
initial  review  of  fans  available  from  the  unit  heater  manufacturer.  New  York  Blower 
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Company,  indicates  that  fans  with  larger  flow  rates  having  similar  dimensions  as  those 
presently  installed  are  available.  An  excerpt  from  a  New  York  Blower  technical 
publication  showing  various  size/flow  rate  fans  is  included  in  Appendix  E. 

2.  Increasing  Entering  Water  Temperature 

a.  Analysis  Method 

The  hydronic  system  analysis  performed  thus  far  was  done  assuming  an 
entering  water  temperature,  T^,  of  190°F,  an  industry  accepted  hydronic  system  operating 
temperature.  It  should  be  noted  that  a  non-pressurized  water  system  will  never  match  the 
heat  transfer  performance  of  a  steam  system  with  all  other  system  parameters  equal.  This 
is  due  to  the  difference  in  the  water’s  and  steam’s  entering  temperature,  which  as  was 
shown  previously,  directly  affects  the  Effectiveness-NTU  computation  of  the  maximum 
heat  transfer  rate,  q  which  in  turn  affects  the  heat  transfer  rate,  A ,  (recall  that 

^max  ^ 

6  ^  C  .  (T-TA  and  Q  =  eQ  )• 

b.  Analysis  Results 

The  effects  of  increasing  the  inlet  water  temperature  from  190°F  to 
205  °F  are  shown  in  Tables  15  and  16. 

c.  Practicality 

The  means  to  provide  an  increased  water  temperature  are  available  from 
many  commercial  sources.  An  engineering  judgement  needs  to  be  made  as  to  whether  the 
energy  put  into  additionally  heating  the  water  is  worth  the  energy  received  in  the  form  of 
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Heater 

^7’3  =  190°F 

(BTU/hr) 

Qt^=205°F 

(BTU/hr) 

^r3=205°F 
^r3  =  190®F 

B-25 

11,713 

13,064 

1.11 

B-70 

28,214 

31,469 

1.11 

Type  B-25  include  unit  heaters  2-8-0  and  2-73-1 . 

Type  B-70  include  unit  heaters  2-80-1, 2-40-1, 2-40-2, 2-60-1, 2-60-2. 


Table  15:  Unit  Heater  Hydronic  Performance  Data  with  Increased  Entering  Water  Temp. 


Heater 

S7’3  =  190°F 

(BTU/hr) 

Ar 

^*''r3  =  190‘’F 

(°F) 

Qtj=2Q5°F 

(BTU/hr) 

At 

=205  “F 

(°F) 

Qt^=20S°F 

At 

^*''r3  =205“F 

^7’3  =  190°F 

At  . 

01-24-1 

74,778 

19 

83,541 

21 

1.11 

1.11 

01-25-1 

66,393 

44 

70,920 

47 

1.07 

1.07 

01-50-0 

43,215 

53 

46,162 

57 

1.07 

1.07 

2-25-1 

22,526 

20 

24,779 

22 

1.10 

1.10 

2-16-1 

15,152 

40 

16,720 

44 

1.10 

1.10 

Table  16:  Duct  Heater  Hydronic  Performance  Data  with  Increased  Entering  Water  Temp. 


heated  air.  A  control  system  that  would  allow  water  temperatures  to  be  easily  adjusted  to 
higher  temperatures  in  colder  weather  is  desirable.  To  minimally  impact  a  hot  water 
heating  system  on  the  whole,  local  electrical  booster  heaters  could  be  installed 
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near  heat  exchanger  water  inlets.  This  could  be  a  desirable  means  to  provide  hotter  water 
to  the  duct  preheaters,  which  directly  heat  weather-intake  air. 

E.  KEY  FACTORS  IN  HEAT  EXCHANGER  PERFORMANCE 

1.  Dependence  on  Air  Side  Convection  Coefficient 

Results  of  Effectiveness-NTU  analysis  for  both  the  hydronic  and  steam  system 
heat  exchangers  are  dependent  on  the  air  side  convection  coefficient,  h„j^.  Retracing  the 
path  in  obtaining  heat  transfer  rates  and  air  temperature  rises  from  the  Effectiveness-NTU 
analysis,  it  was  discerned  that  the  most  influential  characteristic  of  a  given  heat  exchanger 
was  the  air  side  convection  coefficient.  Due  to  the  its  relatively  small  magnitude,  the 
value  of  dominated  computation  of  the  overall  coefficient,  From  was 
determined  the  number  of  transfer  vmits,  NTU.  Larger  values  of  NTU  yielded  larger 
efficiencies,  e,  which  when  multiplied  by  the  maximum  heat  transfer  rate,  ,  yielded 


increased  heat  transfer  rates,  Q . 

Gas-liquid  heat  exchanger  performance  is  thus  regarded  as  air  side  limited,  with 
the  air  side  convection  coefficient  the  key  parameter.  The  manner  in  which  is 
computed  is  worth  revisiting  to  determine  influencing  factors.  Recalling  from  Chapter 
III: 


(37) 
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In  this  expression,  the  use  of  Tj  (inlet  air  temperature,  the  lowest  temperature  in 
the  system)  results  in  the  largest  difference  when  subtracted  from  Tg,  thus  giving  the  most 
conservative  (lowest)  values  of  Carrying  this  conservative  value  of  through  the 
analysis  as  described  above  results  in  conservative  values  for  heat  transfer  performance. 

2.  Manufacturer’s  Stated  Heat  Transfer  Performance 

Using  a  manufacturer’s  air  temperature  rise  table,  such  as  the  on  shown  in  Figure 
2,  the  steam  system  heat  transfer  rate  can  be  solved  for  directly  ( i.e.  Q  =  mc^[T^-T^]). 

Proceeding  to  solve  for  e,  NTU,  and  (neglecting  all  thermal  resistances  but  air 
side  convective  resistance  for  simplification,  i.e.  ^lUA  =  will  give  an 

indication  of  what  T ,  equivalent  a  manufacturer’s  resulting  heat  transfer  rates/air 

temperature  rises  are  based  on. 

Looking  at  manufacturer’s  information  presented  in  Figure  2,  and  considering  the 
case  of  a  single  row,  12  fin  per  inch  heat  exchanger,  with  an  air  flow  velocity  of  500 
FPM,  the  temperature  rise  from  O^F  is  77.2 °F.  Solving  for  and  then  the  T,  equivalent 
as  outlined  above  yields  a  value  of  7)  equal  to  58°F  (see  Appendix  F  for  computation 
details).  This  is  in  contrast  to  the  value  of  0°F  used  for  computations  of  /z^j,  in  the 
analysis  done  in  Chapter  III  (see  output  for  Fortran  computer  code  “hxair”  I  Appendix  B). 
This  demonstrates  that  values  of  /z„,y  that  manufacturer’s  heat  transfer  rates  are  based  on 
are  significantly  higher  than  values  that  were  used  in  this  thesis,  leading  to  differences  in 
heat  exchanger  performance  values. 
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The  difference  in  the  analysis  starting  point  oiT,  outlined  above  indicates  that 
magnitudes  of  heat  exchanger  performance  data  found  in  this  thesis  are  conservative  by 
approximately  15-20  percent  when  compared  to  manufacture’s  data.  The  percent 
difference  between  the  hydronic  and  steam  system  heat  transfer  performance  will 
however  still  be  consistent  with  the  results  of  this  thesis  -  compared  data  will  have  larger 
magnitudes  but  with  approximately  the  same  percent  differences.  The  only  true 
verification  of  any  prediction  of  the  complex  heat  transfer  involved  in  a  compact  gas- 
liquid  heat  exchanger  is  through  controlled  measurements  under  known  conditions.  By 
experimental  verification,  predictions  to  cover  similar  geometries  are  then  able  to  be 
made  with  more  certainty. 

F.  OPTIMIZATION  OF  WATER  FLOW  RATES 

1.  Optimization  Candidates 

Heat  exchangers  that  meet  or  exceed  the  heat  transfer  rates  air  temperature  rises  as 
stated  in  the  ship’s  prints  (heat  exchangers  2-8-0, 2-80-1,  and  01-24-1)  can  be  further 
analyzed  to  determine  an  optimal  water  flow  rate.  With  the  air  flow  rate  in  this  view  of 
system  operation  fixed  at  the  manufacturer’s/print  values,  the  water  mass  flow  rate  can  be 
adjusted  to  lesser  values  to  obtain  an  optimal  value  at  which  the  minimum  heat  delivery 
rate  and  air  temperature  rise  can  be  achieved.  The  objective  in  this  optimization  is  to 
achieve  the  minimum  water  pumping  power  required  to  meet  stated  system  heat  transfer 
performance  parameters. 
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2.  Effectiveness-NTU  Analysis  Method  Readapted 
The  Effectiveness-NTU  analysis  method  can  be  used  once  again  with  the  intention 
of  solving  for  the  heat  exchanger  surface  area  for  a  given  heat  transfer  rate.  The  heat 
exchanger  surface  area  can  be  solved  for  using  a  relation  presented  in  Chapter  III,  i.e. 


(38) 


where,  for  instances  of  (m c  )  ,  <  (me  )  .  : 

^  ^  p'wtr  p'air 


(39) 


Thus,  for  a  given  heat  exchanger  with  known  air  flow  rate,  air  side  convection 
coefficient,  estimated  water  side  coefficient,  and  stated  heat  transfer  rate,  the  air  side 
surface  area  can  be  computed  and  compared  with  the  known  value.  This  optimization 
procedure  is  continued  until  the  computed  surface  area  is  less  than  or  equal  to  the  known 
surface  area. 

Once  an  optimal  water  mass  flow  rate  is  found,  the  accompanying  water  side 
convection  coefficient  is  adjusted  to  match  the  value  used  in  the  computation  of  the 
overall  coefficient.  This  is  not  too  onerous  of  a  task  since  as  was  described  in  the  water 
side  convection  computation  section  of  Chapter  III  (results  shown  in  Appendix  C),  the 
water  side  convection  coefficients  at  low  water  flow  rates  (i.e.  0.5  gpm)  are  already  an 
order  of  magnitude  above  the  air  side  convection  coefficients  and  thus  have  little  effect 
on  the  overall  coefficient. 
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a.  Sequence  of  Computations 

Computations  involving  the  optimizations  of  water  flow  rates  using  the 
Effectiveness-NTU  analysis  were  completed  with  assistance  of  Fortran  computer  codes. 
The  sequence  of  computations  in  the  Fortran  code  titled  “hxNTUOPT.f shown  as  part 
of  Appendix  C,  is  summarized  in  Figure  24.  Program  inputs  are  the  previously 
computed  estimated  values  of  water  inlet  and  outlet  temperatures  Tj  and 
respectively,  air  volumetric  flow  rates,  heat  exchanger  air  side  surface  areas  A  ,  and  the 

^air 

stated  design  heat  rates  n 

Sprints 

b.  Tabular  Results 

An  example  of  Effectiveness-NTU  water  mass  flow  rate  optimization 
program  output  for  heater  2-80-1  is  shown  in  Table  17.  Program  outputs  for  other 
optimized  heat  exchangers  2-8-0  and  01-24-1  are  shown  in  Appendix  C.  The  shaded  last 
row  represents  the  final  optimized  parameters. 

A  summary  of  optimized  water  mass  flow  rates  for  unit  heaters  2-8-0, 2- 
80-1,  and  duct  heater  01-24-1  are  shown  in  Table  18. 

c.  Low  Optimal  Flow  Rates 

As  shown  in  Table  1 8,  for  heat  exchangers  2-8-0  and  2-80-1,  the  minimum 
water  mass  flow  rate  are  very  low  compared  to  the  value  of  10  gpm  used  for  other  heat 
exchangers  in  the  system.  It  is  apparent  that  these  two  heat  exchangers  could  have  much 
less  heat  transfer  surface  area  on  both  the  water  and  tube  sides  at  the  expense  of  an 
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increased  water  mass  flow  rate  and  still  provide  heat  transfer  rates  stated  on  the  prints. 
This  can  be  seen  from  the  results  in  Table  17,  where  the  increase  in  water  mass  flow  rate 
corresponds  to  a  decrease  in  air  side  surface  area.  Water  flow  rates  as  low  as  those 
obtained  for  heat  exchangers  2-8-0  and  2-80-1  delivering  heat  transfer  rates  stated  on  the 
ship’s  prints  indicates  that  the  heat  exchangers  are  oversized  for  the  application.  This 
apparent  oversizing  is  also  evident  when  it  is  noted  that  the  same  size  unit  heater  (B-70) 
is  used  to  deliver  heat  transfer  rates  varying  from  5,294  to  64,860  BTU/hr.  Selections  of 
apparently  oversized  heat  exchangers  for  a  particular  heat  load  was  most  likely  based  on 
other  sound  reasoning,  not  solely  on  heat  transfer  capacity. 


(Ibm/hr) 

V  , 

yiftr 

(gpm) 

C  C 

min  max 

(BTU/hr-R) 

C 

r 

e 

NTU 

Computed 

A 

‘‘air 

(in^) 

59 

0.12 

59.6 

431.6 

0.14 

1.00 

5.00 

10,216 

69 

0.14 

69.6 

431.6 

0.16 

0.86 

2.16 

5,156 

79 

0.16 

79.6 

431.6 

0.18 

0.75 

1.52 

4,150 

89 

0.18 

89.7 

431.6 

0.21 

0.66 

1.20 

3,689 

99 

0.20 

99.7 

431.6 

0.23 

0.60 

1.00 

3,418 

109 

0.22 

109.7 

431.6 

0.25 

0.54 

0.86 

3,235 

119 

0.24 

119.7 

431.6 

0.28 

0.50 

0.75 

3,079 

129 

0,27 

\ 

lilliil 

0,30-  = 

0,4^ 

0,67  : 

Table  17.  Results  of  Water  Mass  Flow  Rate  Optimization  Using  Effectiveness-NTU 
Analysis  for  Heater  2-80-1.  Actual  A  =  3,013  inland  Qprmts  -'I  STUlhr  _ 

^  air 
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Heater 

Number 

Q 

^  prints 

(BTU/hr) 

Actual 

(in^) 

Water  Flow  Rates  at 

A  {computed)  <  A  {actual) 

^air  ^air 

(Ibm/hr)  (gpni) 

2-8-0 

7,747 

3,013 

129 

0.27 

2-80-1 

5,294 

5,198 

51 

0.10 

01-24-1 

51,600 

17,700 

812 

1.67 

Table  18.  Summary  of  Effectiveness-NTU  Optimization  Analysis  Results 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  SYNOPSIS 

This  thesis  successfully  analyzed  heat  transfer  performance  aspects  of  a  set  of  heat 
exchangers  installed  on  the  U.S.  Coast  Guard’s  WTGB  Icebreaking  Tug  class  cutters. 
Initial  analysis  with  acknowledged  conservative  definitions  of  air  side  convection 
coefficients  determined  that  the  hydronic  system  provided  on  average  seventy  percent  of 
the  heat  transfer  capabilities  available  with  the  steam  system.  Practical  improvements  to 
the  hydronic  system  were  shown  to  increase  heat  exchanger  performance  parameters  by 
an  average  of  ten  percent.  It  was  notable  that  the  added  heat  transfer  available  from  steam 
is  not  due  to  a  property  of  steam  itself  such  as  latent  phase  change  effects,  but  is  due 
solely  to  the  increase  in  entering  tube  side  temperature.  Judging  by  heat  transfer 
capabilities  alone,  with  the  described  conservative  assumptions  on  which  these  results  are 
based,  use  of  currently  installed  heat  exchangers  in  a  hydronic  system  is  a  viable  option. 

B.  PROPOSED  HYDRONIC  SYSTEM 

A  hydronic  system  circulating  water  at  10  gallons  per  minute  and  190  degrees  F, 
will  provide  at  minimum  the  heat  transfer  and  air  temperature  rises  shown  in  Table  10. 
These  performance  parameters  were  arrived  at  considering  the  entire  thermal  circuit 
involved  and  with  conservative  assumptions  (when  compared  to  manufacturer’s  methods) 
for  air  side  convective  resistance  determinations.  Hydronic  system  performance 
parameters  can  be  improved  through  the  practical  methods  outlined.  Considering  heat 
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transfer  capabilities,  use  of  the  currently  installed  heat  exchangers  as  part  of  a  hydronic 
system  is  plausible. 

Use  of  a  pressurized  water  system,  where  system  pressure  would  allow  the 
circulating  water  temperature  to  equal  or  exceed  that  of  steam  (240  °F)  is  another 
approach  to  the  hydronic  system  design  that  could  be  explored. 

C.  FURTHER  STUDY 

There  are  numerous  issues  to  examine  further  in  order  to  take  into  account  all 
aspects  of  using  presently  installed  heat  exchangers  as  part  of  a  hydronic  system.  Issues 
remaining  are  related  to  the  hardware  involved  as  well  as  further  analysis/design  work. 

One  hardware  issue  remaining  is  the  internal  heat  exchanger  component’s 
compatibilities  with  circulating  hot  water  as  opposed  to  steam.  In  steam  heat  exchangers, 
there  is  typically  a  strip  of  shaped  metal  which  distributes  the  flow,  insuring  that  the  fluid 
passes  nearly  equally  through  all  tubes.  The  erosion  effects  of  water  on  this  and  other 
internal  elements  designed  for  use  with  steam  should  be  examined. 

A  hardware  issue  on  the  external  side  of  the  heat  exchanger  is  the  relative  location 
of  water  vs.  steam  inlet  and  outlet  piping.  One  manufacturer’s  information  indicated  that 
for  a  steam  heat  exchanger,  inlet  piping  is  above  the  outlet  to  allow  for  condensate 
drainage.  Conversely  for  circulating  water,  the  inlet  is  below  the  outlet  to  help  prevent  air 
in  the  system  from  inhibiting  water  flow.  The  effort  involved  in  possibly  reversing  inlet 
and  outlet  piping  when  switching  from  a  steam  to  hydronic  system  thus  needs  to  be 
examined. 
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Further  design  and  analysis  work  includes  among  other  matters,  the  hydronic 
system  design  to  provide  water  at  either  10  gallons  per  minute  or  the  optimized  values  to 
the  respective  heat  exchangers.  Pressure  drops  across  each  heat  exchanger  and  across 
other  system  components  need  to  be  accurately  gauged  to  arrive  at  a  required  pumping 
head.  Air  side  pressure  drops  also  need  to  be  examined  if  any  of  the  heat  transfer 
enhancements  involving  increases  in  air  flow  rate  are  adopted,  to  insure  fans  remain 
properly  sized. 
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APPENDIX  A.  WTGB  CLASS  CUTTER 
ILLUSTRATION 
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APPENDIX  B.  HEAT  EXCHANGER  DIMENSIONS 
AND  MATERIAL  PROPERTIES 
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New  York  Blower 

Company® 

7660  QUINCY  STREET— WILLOWBROOK,  ILLINOIS  60521 


AVb‘ 


SUPPLY  AND  RETURN 
FEMALE  PIPE  TAP 


UNCASED  COILS 


DIMENSIONS  (INCHES) 


Size 

Coil 

Face 

Area 

A 

B 

D 

2 

.70 

13 

10 

10 

46 

.96 

lays 

llVz 

11% 

812 

1.29 

ISVs 

13V2 

13% 

1420 

1.88 

19% 

iev2 

16% 

1824 

2.78 

22% 

20 

20 

4256 

4.92 

31 

27 

26% 

y  y 


DIMENSIONS  [INCHES] 


Size 

A  or  B 

A 

B 

C 

max. 

D 

[FPT] 

Wheel 

dia. 

Approx.* 

weight 

[lbs.] 

25 

12 

13‘/4 

2078 

143/4 

172 

8 

65 

45 

13% 

ISVa 

2078 

1578 

172 

10 

75 

70 

153/4 

1678 

2072 

17% 

172 

12 

115 

105 

I8Y4 

1972 

2072 

21 

2 

14 

145 

120 

135 

155 

22 

2378 

2474 

2478  ' 

2 

18 

180 

180 

200 

200 

240 

270 

300 

2874 

3174 

243/8 

323/4 

272 

24 

305 

305 

305 

310 

"Weights  will  vary  with  motor  specifications.  Tolerance:  ^Vs" 
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COLMAC 


ICOIL 


Manufacturing  Inc. 


STEAM  HEATING  COILS 
5/8”  OD  TUBE 


Dimensions  For  Steam  Coils 


5/8"  Basic,  Heavy  Doty  Steam  Heating  Coil 


NDTEt 

1.  5/16'  DIA  HOLES  ON  3'  CTRS  FROM  CENTERLINE  OF  CASING. 

a  COIL  PITCHED  IN  CASING'  TOWARD  RETURN  END  1/4'  PER  FOOT  OF  FINNED  LENGTH. 
B  DIM  PLUS  C  DIM  PLUS  FIN  HEIGHT  EQUALS  H. 


88 


Carrier 


UNiT  SIZE  39B 


aSMBMI 

400 

500 

600 

700^ 

Foce  Velocity  (^m) 

800 

(H  &  V  units  only) 

900 

1000 

1100 

_ 

1200 

040 


Fqcc  Area  (sq  ff) 

U>B  end 

Steam  Distributing  Tube 

“humidifier  conn: - 

(no.  ...size) 

Atom!  zing  Spray 
Steom  Grid 
Supply  (in.  OD) 

Drain  (in.  OD) 

RATING  WT  (lb)t - - 


(  UJU 

V 

861 

1255 

1076 

1568 

1291 

1882 

1507 

2196 

1720 

.2512 

1935 

2826 

2150 

3140 

2365 

3454 

i . .  2580 

3768 

2.15 

1 

3.14  1 

1.73 

2.72  ; 

060 

2142 

2677 

3^13 

3749 

4288 

4824 

5360 

5896 

6432 


5.36 

4.81 


FACE  AND  BYPASS 

damper  areas  (sq  ft) 

MIXING  BOX 
PAMPER  AREAS  (sq  ft) 


450 


6.0 


7.1 


PANS  (no.  '*4ic«l  I., .no.  outlets) 
|Diom**(in.) 
Length  (in.) 

I  No.  of  Blodes 
Shoft  Critical  Speed  (rpm) 

Max  Operating  Speed  (rpm) 
Fan  Sheave  Bore  (in.) 

Max  Broke  Hp 


5.9 


U.O 


6.9 


1...1 

9y, 

7y, 

43 

2700 

2160 

IX. 


10.4 


1...1 

lOV, 

8  ■ 
48 

2550 

2040 

IX. 


1...I 

i2y, 

9S 

43 

1700 

1360 


base  unit 

fan  &  COIL  section  l 

w 

_  H 

Mi;tlNG  BOX  SECtlOTi - T” 

Type  A 

46y, 

38^ 

34^ 

— 

DSSEI 

m 

HiiCiUl: 

tSHig 

H 

tSBi 

BOS 

tWill.l 

m 

nsQ2Q]i 

mg 

B 

Pm 

|H  &  V 
33V 
42y, 
23V 

_ 

fType  A 

67V 

46V 

43 

c 

jType  B 
49V 
46V 

i« 

Type  D 

28  V 
46V 

X 

41  4 
36*/, 
25V, 

— dV — 

dir, 

36’/, 

25V 

- aU 

21V 

36V 

25V 

1  21V 

36V 

lj5V, 

24V 

36V 

29V 

^4V 

36*/, 

29V 

“T4V 

36V 

29V 

’14^ 

36V 

29V 

'30V 

44V 

37V 

30V 

44V 

37V 

70 
30V 
44  V 

1  4vr^ 

30V 

44V 

FILTER  SECTION  W 

H 

preheat  SECTION - [T 

36?; 

20% 

4^ 

36*/, 

20V 

4V 

36V 

20V 

36V 

20V 

4V 

36*/, 

23V 

4V 

36*/, 

23V 

4V 

36V 

23  V 

36V 

23y 

■B 

4V 

44  V 

J/  /, 

4V 

44  V 

703/ 

J  /  .'B 

4X 

44V 

W 

__  H 

^WasTplInDm"^ - iT 

6 

36V, 

20V 

6 

36V 

20V 

19 

ES 

H 

1^1 

19 

msM 

H 

9 

44V 

29V 

;  zV/, 

I  9 

44V 

29y 

dV  4 

9 

44  V 

703/ 

29V 

FACE  AND  BYPASS  W 

damper.  SECTION  h 

"service  AREA  -^"" - A“ 

3^/, 

28V, 

6 

36V 

28V 

0 

36V 

28V 

- 

M 

WsM 

H 

m 

msm 

■ 

9 

44  V 

A}V 

9 

4dV 

A^r/ 

zv  4  • 

9  1 

44V  1 

I 

- 

requirements  B 

L  -  Length  ^ 

W  -  WirtTh  . 

46V 

48V 

-  H< 

34 

48V 

»tghT 

19V 

48V 

40V 

48V 

13V 

48V 

39V,  ~22% 

^8%  48H 

m 

fSM 

1^^ 

m 

-17 — 
56V 
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APPENDIX  C.  FORTRAN  COMPUTER  CODES 

AND  OUTPUTS 
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APPENDIX  D.  EFFECTIVENESS  -  NTU  ANALYSIS 

TABULAR  RESULTS 
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The 

New  York  Blower 

Company® 

7660  QUINCY  STREET— WILLOWBROOK,  ILLINOIS  60521 


Model  N 


DIRECT-DRIVE  PROPELLER  FANS 

EXHAUST  OR  SUPPLY 

•  Eight  wheel  diameters — 8”  through  24". 

•  243  to  6420  CFM —  up  to  V2*'  static  pressure. 

•  Panels— square  steel  construction  with  streamlined  venturi 
inlet... venturi  is  reversed  in  supply-fan  panels... baked-ereen 
alkyd  finish. 

•  Wheels — aluminum  blades  with  steel  hubs. 

•  Motor  mounts—wi re-guard-type  motor  mount  [see  photo  at  left] 
is  standard  on  all  direct-drive  units... guard  is  zinc-plated  steel. 

•  Motors  standard  motors  are  totally  enclosed  air  over  with  pre- 
lubricated  ball  bearings  except  1/12  and  V20  HP  motors,  which 
are  shaded-pole  totally  enclosed  permanently  lubricated  sleeve- 
bearing  type. 

Motors  1/4  HP  and  larger  are  suitable  for  either  horizontal  or 
vertical  service... specify  “for  vertical  mounting”  to  have  wheel 
locked  to  motor  shaft. ..  1/20  and  V12  HP  motors  are  not  suitable 
for  vertical  service. 


H  =  min.  wall 
thickness. 


SPECIFICATIONS  DIMENSIONS  IN  INCHES. 


EXHAUST 


NOTEf&hSntoTillSe'J^  ""0  PacKaging. 

Whenordenng.spec,fyconnp.etemode,numbe,asshownonpage3.  D^Se^nsions  o'?,  Tolaranc, 


H  minimum 

Mounting  hole 

Et 

Gt 

Auto* 

Motor 

no.  and  diameter 

Weight^ 

[lbs.] 

101/4 

101/4 

15/8 

41/2 

8-5/16 

8-9/32 

25 

103/8 

101/8 

15/8 

43/4 

8-5/16 

8-9/32 

29 

111/2 

103/4 

2 

53/8 

8-5/16 

8  -  9/32 

35 

111/4 

lOVe 

2 

51/8 

8-5/16 

8-9/32 

40 

12 

11 

2 

51/2 

8-5/16 

8-9/32 

50 

111/2 

107/8 

2 

51/8 

8-5/16 

8-9/32 

65 

127/8 

12 

2 

53/8 

8-5/16 

8  -  9/32 

80 

133/8 

123/8 

2 

51/2 

8-5/16 

8-9/32 

95 

101/4 

101/4 

91/2 

8-5/16 

8-9/32 

25 

101/4 

101/4 

Auto¬ 

91/2 

8-5/16 

8  -  9/32 

29 

1 1 

11 

matic 

1  9V2 

8-5/16 

8  -  9/32 

35 

11  Vs 

11V8 

shutter 

91/2 

8-5/16 

8-9/32 

40 

111/4 

111/8 

not 

91/2 

8-5/16 

8  -  9/32 

50 

12 

12 

available 

91/2 

8-5/16 

8  -  9/32 

65 

121/2 

121/2 

121/2 

121/4 

91/2 

91/2 

8  -  5/16 
8-5/16 

8  -  9/32 

8-9/32 

80 

95 

Tolerance;  ±  Ve’ 
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APPENDIX  F.  COMPUTATION  OF  AIR  SIDE 
CONVECTION  COEFFICIENT  USING 
MANUFACTURER’S  DATA 
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